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Abstract 

A simple generalization of the MHD model accounting for the fluc- 
tuations of the configurations due to kinetic effects in plasmas in short 
times small scales is considered. The velocity of conductive fluid and the 
magnetic field are considerd as the stochastic fields (or random trial tra- 
jectories) for which the classical MHD equations play the role of the mean 
field equations in the spirit of stochastic mechanics of E. Nelson. 



PACS codes: 47.10ad, 47.27.E-, 47.27.ef 
Keywords: 



1 Introduction 

Theoretical investigations of cross-field transports in the operating ITER- 
FEAT (International Thermonuclear Experimental Reactor) calls for an increas- 
ing confidence in the modelling efforts that force one to search for the new prin- 
ciples of simulations. The aim of our work is to provide a possible ground for 
the optimization of existing numerical simulation algorithms for the large scale 
simulations in hydrodynamics. 

The set of equations which describe magnetohydrodynamics (MHD) are a 
combination of the Navier-Stokes equations (jTJ 1) of fluid dynamics and Maxwell's 
equations (jTj2) of electromagnetism, 

v + VP + (v ■ V)v — (b ■ V)b = i/Av 
b+{v- V)b - (b ■ V)v = AAb, ' ! 
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in which A = c 2 /4ira is the resistivity constant, the inverse Prandtl-type con- 
stant. The equations ([1]) describe the dynamics of electrically conducting fluids. 
The fluid velocity v(x,t), x S M. d is supposed to be incompressible, so that 
the mass continuity equation for that is reduced to the transversally condition 
V • v = 0. The normalized magnetic field h(x, t) is 

b(x,t) = ^d, V-b = 0, (2) 

in which B(x,i) is the magnetic induction and p is the density of medium. 
In addition to the usual hydrodynamical interaction presented in the Navier- 
Stokes equation, in the moving medium, there is the Lorentz force pQ exerted 
on charged particles in the electromagnetic field, 

[curl B x B] = (B ■ V)B - V(B 2 /2). (3) 

The first term in ([3]) is amended to the hydrodynamical interaction while the 
second term redefines the pressure field in the medium, 

B 2 

P(x,t)^p(x,t) + —. (4) 

The dynamical equation describing the evolution of magnetic field follows from 
the simplest form of Ohm's law pQ, 

i(x, t)=a (e(x, t) + ~ [<p(x, t) x B(x, t)] , \ (5) 

in which a is the conductivity, c is the speed of light, E and tp are the electric 
field and the electrostatic potential respectively, and the Maxwell equations 
neglecting the displacement current. 

The purely longitudinal contributions of the pressure gradient VP and of 
interactions can be eliminated from ([!]) by the applying of transverse projection, 

if written in the Fourier space. 

The equations (p} are applicable to the plasma if it is strongly collisional, 
so that the time scale of collisions is shorter than the other characteristic times 
in the system, and the particle distributions are therefore close to Maxwellian. 
When this is not the case (for instance in fusion plasmas), we are interested in 
smaller spatial scales, in which it may be necessary to use a kinetic model which 
properly accounts for the non-Maxwellian shape of the distribution function. 
However, because MHD is very simple, and captures many of the important 
properties of plasma dynamics, it is often qualitatively accurate, and therefore 
while accounting for the possible kinetic effects in plasmas we are nevertheless 
interested to stay within the general framework of the MHD approach. 
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In the present paper, we investigate a simple generalization of the MHD 
model (Q]) modelling fluctuations of the configurations {v, b} due to kinetic ef- 
fects in plasmas. In the proposed model, we suppose that {v, b} are the stochas- 
tic fields for which |T]) plays the role of the mean field equations. Recently, we 
have implemented the similar approach for the Burgers and Navier-Stocks equa- 
tions in 



2 Stochastic dynamics as the Brownian motion 

It is well known that many problems in stochastic dynamics can be treated as 
a generalized Brownian motion (6 (u — u (x, t)))^, in which the classical random 
field indicating the position of a particle u(x, t) meets a Langevin equation, 

u(s,t) = K(u) + Q(u)+£, (7) 

where £ is the Gaussian distributed stochastic force characterized by the corre- 
lation function 

D e = (£0- 

Here the angular brackets (. . .), denote an average position of particle with 
respect to the statistics of £. K(u) is the linear differential operator, and Q(u) 
is some i-local (independent of time derivatives) nonlinear term which depends 
on the position u(x, t) and its spatial derivatives. Such a representation was a 
key idea of the famous Martin-Siggia-Rose (MSR) formalism, [3]-[6]. 

An elegant way to obtain the field theory representation of stochastic dy- 
namics is given by the functional integral 

(«(«-u(x,t))) c = J D{[exptr(~l;DAs(u-u(x,t)) (8) 

where the tr-operation means the integration J dxdt and the summation over 
the discrete indices. The instantaneous positions u (x, t) meet the dynamical 
equation that can be taken into account by the change of variables 

<J(u-u(x,*))-K$(u(x,t)-Q(u)-£) (9) 

should the solution of dynamic equation exists and is unique. The use of integral 
representation for the 8— function in ([8]) transforms it into 

J D^DuDu' exp tr f-^D^-u'u + u'Q(u) +u'A detM, (10) 

in which u' (x, t) is the auxiliary field that is not inherent to the original model, 
but appears since we treat its dynamics as a Brownian motion. The Jacobian 
detM relevant to the change of variables ([9]) is discussed later. 

The Gaussian functional integral with respect to the stochastic force £ in 
(fT0|) is calculated 

J DuDu' expS (u,u r ) detM, (11) 
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in which 

S(u,u') = tr 



—u'D^u' — uii + u'Q (u) 



(12) 



By means of that all configurations of £ compatible with the statistics are taken 
into account. The integral (TTUJ) identifies the statistical averages (. . .)^ with the 
functional averages of weight exp S. The formal convergence requires the field 
u to be real and the field u' to be purely imaginary. 

The functional averages in (fTTj) can be represented by the standard Feynman 
diagram series exactly matching (diagram by diagram) the usual diagram series 
found by the direct iterations of the Langevin equation averaged with respect to 
the random force This fact justifies the use of functional integrals in stochastic 
dynamics at least as a convenient language for the proper diagram expansions. 

The Jacobian det M in l|10p depends upon the nonlinearity Q (u) . If Q (u) 
does not depend upon the time derivatives, all diagrams for det M are the cycles 
of retarded lines A cx 6 (t — t') and equal to zero excepting for the very first 
term, 

detM = const ■ exp trA, (13) 

Then, the convention is used for the Heaviside function of zero argument, 9(0) = 
0, so that detM = const, [7]. 

The functional averages computed with respect to the statistical weight 
exp S can be expanded into the series of Feynman diagrams drawn with the 
interaction vertices determined by the nonlinearity Q(u) and two propagators 
(lines) which have the following analytical representations (in the Fourier space) 

A uu ' — "7 ^"\~> A uu — "/ I \~; (14) 



-iuj + K 



(uj 2 + k 2 ^ 



in which K is the Fourier image of the linear part in the Eq.©. The inverse 
Fourier transform of A uu > shows that it is retarded, A uu i cx 6(t — to). 

Many dynamical systems are driven by the non-random external forces. It is 
worth to mention that if one assumes |£| — > (and consequently — > 0), then 
the Feynman diagram series is trivial since the propagator vanishes, A UM = 0. 
However, the diagram series would be recovered by means of regular external 
forcing [2\ that gives rise to a branching representation of stochastic dynamics. 

3 Probabilistic interpretations for the solutions 
of elliptic equations 

The striking similarity between the Schrodinger equation written for free 
particles and the diffusion equation motivated the search for a stochastic in- 
terpretation of the quantum mechanics. The first attempt had been made by 
E. Schrodinger himself [8] and accomplished by J.C. Zambrini who derived the 
genuine Euclidean version of quantum mechanics [9 . 
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As a counter motion, E. Nelson [TO] had proposed a generalization of the 
theoretical scheme of classical mechanics known as stochastic mechanics. Clas- 
sical deterministic trajectories are substituted by random trajectories of well 
defined stochastic processes. Under appropriate conditions, and for a large 
class of dynamical systems, the basic equations of stochastic mechanics show a 
surprising connection with the basic equations of quantum mechanics, |llj-|12j. 
Therefore, stochastic mechanics gives an approach to quantization of dynami- 
cal systems, based on methods of probability theory and stochastic processes. 
The original formulation of stochastic mechanics rests on two basic hypothesis, 
[13]- [2]. The first assumes that the trajectories of the dynamical system are 
perturbed by an underlying Brownian motion. The second is a particular form 
of the second principle of dynamics, where the classical acceleration is replaced 
by a suitable form of stochastic acceleration. Further developments of the the- 
ory show that the basic equation of stochastic mechanics can be derived from 
variational principles, in complete analogy with classical mechanics, based on 
the same classical action, but exploiting stochastically perturbed trajectories as 
trial trajectories [JJ]. The basic equations of stochastic mechanics (the conti- 
nuity equation and Madelung equation) can be immediately connected with the 
Schroedinger equation of quantum mechanics. The entire operator structure 
of the quantum mechanical observables can be easily derived from the general 
structure of stochastic mechanics. From this point of view, stochastic mechanics 
can be considered as a kind of probabilistic simulation of quantum mechanics, 

Stochastic mechanics can be based on variational principles of Lagrangian 
type. In order to obtain that it is necessary to generalize the action of classic 
mechanics to the case where the trial trajectories belong to stochastic processes. 
This program has been partially realized in [T5]-[IE] for dynamical systems on 
curved manifolds. This was the analog of the well known problem of writing the 
Feynman path integral for a quantum system on a curved manifold, so that all 
the results could be immediately translated into the language of Feynman path 
integrals. 

Independently of quantum mechanics and quantum field theory, a proba- 
bilistic interpretation for the solutions of linear elliptic and parabolic equations 
with Cauchy and Dirichlet boundary conditions had been proposed [17] • A 
stochastic process had been defined for which the mean values of some func- 
tionals coincide with the solution of the deterministic equations. The problem 
of existence of such the probabilistic representations for the certain classes of 
nonlinear equations has been studied extensively by Dynkin [18] . The proba- 
bilistic representations of the Fourier transformed Navier-Stokes and Burger's 
equations had been discussed in [TpJ and later extensively developed in the works 
of Oregon group [30]- (S2]. A stochastic representation for the Poisson-Vlasov 
equation has been derived in [33] recently. In all cases when the appropriate 
stochastic process has been constructed, its mean values is the solution of the 
mean field equation, but the process itself always contains more information 
then that of physical relevance in particular. 

A model in which the classical deterministic trajectories u(x, t), x £ M. d , sat- 
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isfying the hydrodynamics equations are substituted by the random trajectories 
of a generalized Brownian motion over the space of fluid velocity configurations 
u driven by the stochastic force £ is known as stochastic hydrodynamics [23] . 
Here, £ is the Gaussian distributed stochastic force characterized by the corre- 
lation function — (££), and the angular brackets (. . .) ? denote an average 
velocity of particle with respect to the statistics of £. In such a formulation, 
the above problem is equivalent to that one of Nelson's stochastic mechanics, 
, |13j . The relevant variational principle leads to an action functional of the 
Martin-Siggia-Rose (MSR) type [3J. The MSR-theory of stochastic hydrody- 
namics has been formulated independently by many authors, [25j-|27j. Diagram 
representations for the Green functions of stochastic hydrodynamics exactly re- 
produce the hydrodynamical diagrams discussed by Wyld, [28]. In general, the 
Green functions of stochastic hydrodynamics diverge for very large moments 
and therefore require the ultraviolet renormalization that has been discussed in 
details in [7j. 

In our previous paper [5] , we have pointed that the diagram representations 
for the Green's functions in hydrodynamics is still nontrivial if one considers a 
regular external forcing, instead of random one. In particular, we have studied 
the S(x — xo)6(t — to) external forcing corresponding to the Cauchy problem of 
the Navier-Stokes equation supplied with an integrable initial condition. Each 
Feynman graph in the diagram series equals to an average over a forest of 
multiplicative branching binary trees (of the certain topological structure) im- 
plemented in [20]- [22]. The branching representations for the Green function of 
Cauchy problem establishes the direct relation between Nelson's stochastic me- 
chanics jllj and the probabilistic interpretations for the solutions of nonlinear 
equations with Cauchy and Dirichlet boundary conditions studied in [18]- [22]. 
It is important to note that in contrast to the MSR theory [7J, the diagrams of 
branching representations for the hydrodynamics equations [2] do not diverge, 
but consitute a regular expansion starting from the standard diffusion kernel. 
Diagram contributions represent the consequent bifurcations of media resulting 
in the cascade of consequent partitions of moments, k = q + (k — q). The 
magnitude of relevant corrections to the standard diffusion spectrum tends to 
zero as t — * to, and the saddle-points (instanton) analysis can be then applied 
to study the "large order" asymptotic contributions [2]. The calculations have 
shown that the asymptotic coefficients demonstrate the factorial growth like 
for the most of models in quantum field theory. The asymptotic series for the 
Green function can be summarized by means of the Borel procedure. In the 
limit t — > to, the corrections to the diffusion kernel have the closed analytical 
form. 
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4 Cauchy problem for the MHD equations and 
its stochastic mechanics formulation 



The Cauchy problem for the MHD equations, 



Vi + 9j (viVj — hbj) — vd 2 Vi = 
bi + dj (vjbi — bjVi) — \d 2 bi = 



5{t - t Q )6{x - x Q ), 
5(t-t' )S(x-x' ), 



(15) 



is supplied with the localized integrable initial conditions vo = v(cc,0) and 
bo = b(x,0). It is possible to construct a stochastic counterpart of MHD (|15[) 
by adding a Gaussian distributed random force into it. This approach leads to 
the stochastic magnetohydrodynamics [59] , [7] , [3D] developed in order to study 
the inertial range scaling laws and different regimes of large-scale asymptotic 
behavior. Alternatively, in the framework of stochastic interpretation, we con- 
sider a stochastic model in which the nonlinear dynamical equations (|15p play 
the role of mean field equations, so that their solutions satisfying the given 
initial conditions play the role of the observables. Instead of the classical deter- 
ministic fields h(x, t) and v(x, t), we study their stochastic trial analogs, b(x, t) 
and v(x,t). Then the Green functions G^, <fi = {v, b} of the original Cauchy 
problem Q150 can be represented by the functional averages, 



G0(x,<;xo,x o ;to,i o ) 



J m expS ($) 



(16) 



over all possible configurations v(x, t) and b(x,t) such that their expectation 
values satisfy (|15[) with the given initial conditions. Here, $ are the functional 
arguments of the action functional such that (jT5|) are its saddle-point 

equations. Its "quadratic" part, Sq(&), corresponds to the linearized equations 
of MHD which play the role of an interaction free theory. We discuss the measure 
of functional integration 2>3> after we consider S(&) . If we introduce the auxiliary 
fields v'(x, t) and b'(x, t), then, up to an inessential constant factor, the relevant 
action functional reads as following: 



S{v,v',b,b') = v'(x ,to) + b'{x' ,t' ) 



-tr 
-tr 



v' i v l + vv' i dj(giv l Vj - gzhbj) - vv'iAv. 
b'ibi + Xb'idj(g3Vjbi - gJ)jVi) - A&'jA&j 



(17) 



where as in (|12p the tr-operator means the integration J dx dt and summation 
over the discrete indices. The quadratic part of (Q~7|) has the form 



Sq(v, v', b, b') = — tr v'iVi + b'ib 



iAvi - XViAbi 



(18) 



The action functionals with ultra- local terms like those presented in (|17|) . v'(xq, to) 
and b'(xa, to), had been studied in |31j . In order to obtain the formal expansion 
parameters in the perturbation theory for the above action functional, we have 
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inserted four coupling constants, 31,2^ = 1, and 33, 4A = 1, in front of the relevant 
interaction terms in (|17|) . Despite their physical dimensions are [51,2] = — M 
and [534] = —[A], their convectional dimensions would be different. If we put 
formally 31,2,3.4 = in (|17p. the action functional turns into S'o( < l > ), free 

of interactions. The correct normalization of integral in (|16p requires that the 
measure of functional integration 23$ be normalized to the volume of orbits in 
the free theory (pure diffusion processes), 

m = ixmimm^ z ^J me ^ So{ ^ (19) 

Moreover, it is obvious that the results of functional averages (TTB1) do not change 
along the set of orbits in the configuration space related by any symmetry trans- 
formation of the action Q17j). The functional integral (|16[) itself is proportional to 
the volume of such orbits. The functional (fP7|) possess the Galilean invariance: 

v(x,<) = v(x + s(t),t) - u(i), 

v'(x,t) = y'(x + s(t),i), 

b(x,t) = b(x + s(i),i), 1 j 

b'(x,i) = b'(x + s(t),t), 

where an integrable function u(t) is a parameter of transformation (the velocity 
of the frame of reference), and its integral s(t) = f^ u(t')dt'. 

The auxiliary fields introduced in (|17|) were not inherent to the original 
physical model, but appear since we treat its dynamics as a Brownian motion. 
While the first two saddle point equations, 

SS „ ss „ ,„ . 

-^ = 0, -^ = 0, 21 

Sv' Sb' 

recover the original Cauchy problem (|15[) (in case 31,2^ = 1, 33, 4A = 1), another 
pair, 

SS SS 

— = 0, = 0, 22 

Sv Sb 

describes the dynamics of auxiliary fields: 

I v'i + dj(givv'iVj - giXb'ibj) = -vAv'i, ^3) 
[ 6'j + dj(g 3 X¥iVj - giw'ibj) = -AA6V 

The above equations are characterized by the negative dissipations since v, A > 
0. Therefore, physically relevant solutions have to satisfy v'(t > 0) = and 
h'(t > 0). The equations (|15I23[) give the time evolution of the infinitesimal 
characteristics of the diffusion processes in MHD and therefore determines them. 
It is important to mention the striking similarity between the equations (|15[ I23|) 
and the coupled system of nonlinear equations for the osmotic velocity and the 
current velocity of stochastically driven Brownian motion, |32j-[35j. 
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We conclude the current section with a remark that the action functional 
(fT7|) can be derived form the standard action functional of the MSR-type for 
stochastic MHD [21], [7] with the ultra- local interaction terms v'(xo,to) and 
h'(x' , t' Q ) added, if one put the stochastic forcing to zero, |£| — > 0. 

5 Diagram technique for stochastic MHD. The 
absence of ultraviolet divergences 

The formal functional averages (|16p computed with respect to the statistical 
weights expS'($) arc interpreted as the infinite diagram series, 

{k 1 +k 2 +k3 + k4 = oc} 

^(a;,t;xo,to;a;o»*o) = fc !fc , fc , fc , 2^ G )> 9i 9 2 9a 9i , 

l ' 2 ' 3 ' 4 ' {*l+*2+AS3+*4 = 0} 

(24) 

where the sum is taken over all positive integer solutions of the equation 

fci + fc 2 + fc 3 + fc 4 = N, h, N e Z+. (25) 

The series (|24p starts from the standard diffusion kernel in the coupling con- 
stants representing the possible physical interactions in MHD: the hydrodynamic 
dragging (51), the Lorenz force (52), the convection (g 3 ), and the stretching 
(34). In view of that, the main tasks of the stochastic approach to MHD is the 
computation of coefficients q^ 1 ^ 2 -^ 3 ^^ j anc [ the estimation of the asymptotic 



properties of the diagram expansions (|24|) . 

In general, the computation {Q < £ 1 > k2 < k3 < ki ) m fl^gj) is a very difficult problem. 
However, if the functional averages fl^|) are taken over the Gaussian distributed 
fields, one can apply Wick's theorem [5S]-[37] to facilitate calculations. In ac- 
cordance to Wick's theorem, the functional averages (HHJ) equal to the sums of 
all complete systems of "pairings" of the interaction operators. If we denote 
operators as nodes and the pairings between them as edges, then each system 
of "pairings" representing the particular average can be visualized by a graph 
called a Feynman diagram. The use of Feynman diagrams helps to make the 
computations visual and clarifies their similarity with the problems of symbolic 
dynamics 38J . 

Any Feynman diagram for the stochastic representation of MHD contains 
one of two ultralocal interaction terms (see Fig. |T|) , with any number of either 
b'— or v'— tails located at the points (x' 0) t' ) and (x ,t ) consequently. They 
reveal the dependence of the responde functions G v and G& in MHD upon the 
initial conditions for the velocity field at (xo,to) and for the magnetic field 
at (x' ,t' ). In quantum field theory, the vertices shown on (fT]) are called the 
composite operators, 

A diagram consists of the " interaction vertices" (see Fig. [5]) , with the factors 
V4s"i Vims' an d Vims i which can be written in Fourier space as 

v v'vv _ 9l v 1 ■ x u , -r 1, \ v v'bb _ 92 V l-s.. u i ,-jf. L. \ 
v ims — 2 \ zm^s i Lu is nj m ) ^ * ims — ^ \ lu im">s t '"js^m ) > 
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Figure 1: The ultra- local interaction vertices in stochastic representation of 
magnetohydrodynamics . 




Figure 2: The interaction vertices in stochastic representation of magnetohy- 
drodynamics. The vertex factors ^ in the Fourier representation are given 
in fl5j5J|. 



K" m b s v = KmS ls k m - ig 3 S im k s ). (26) 

The b'— or v'-tails in interaction vertices (see Fig. [2]) are to be connected with 
the b— or v-tails by means of edges (the response functions of the linearized 
MHD equations): 



Pi(fc) u, Pi(fe) 

ij J- vk 2 ' ij ~~ -iw + Afc 2 ' y ' 



in which the transverse projector P^(fe) is given in The inverse Fourier 
transform with respect to the frequency to in (|27p reveals that the response 
functions of linearized problem are retarded, 

flg" = i/(t-i )P i L(fc)e^ 2 ( t - to ), i^ fc = ff(i-t )P^(fc)e- Afe2 ( t -*°), (28) 

where i? (t) is the Heaviside function supplied by the convention H (0) = 0. It is 
required by the the casualty principle that the time arguments correspondent to 
auxiliary fields v' and b' always precedes the time arguments of v and b. Let us 
note that the inverse Fourier transforms of response functions (|28p correspond 
to the standard diffusion kernels, (in <i-dimensional space), 

*n) 2 



g 4i/(t-to) 

A(g-8 ,*-*o)= , (29) 

V /(4 7 rKi-to)) d 

Similarly to quantum electrodynamics where each Feynman diagram represents 
a certain process of physical interactions between elementary particles, any 
graph in stochastic mechanics drawn with the interaction vertices Figs. [TJ2] con- 
nected by the edges (|2~T|) corresponds to a certain process of interactions between 
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different magnetic and hydrodynamic modes coming along with the cascades of 
consequent partitions of moments, k = q + (k — q). The first diagrams relevant 
to the inverse Green functions (exact response functions) of MHD are shown 
in Fig. [3] and Fig. [?] consequently. Each diagram in expansions corresponds 




Figure 3: The first diagrams for the inverse exact responde function of the fluid 
velocity in MHD. The slash marks auxiliary fields. The first diagram corresponds to 
the responde function of the linearized problem. The second diagram expresses the 
first order correction to the responde function due to the hydrodynamic drag. The 
third diagram stands for the second order correction due to the hydrodynamic drag. 
The forth diagram stands for the first order correction due to the Lorenz force. The 
forthcoming diagrams are related to the corrections risen by the combined effect of 
hydrodynamic and magnetic field fluctuations. 

to a certain analytical integral expression quantifying the contribution of the 
relevant interactions into the entire behavior (they are given in the Appendix 
A). In general case, all integrals with an odd number of factors fej equals zero 
as a consequence of space isotropy. 

The diagram expansion for the Green function would have a definite physical 
meaning if it converges. The standard analysis of ultraviolet divergences of 
graphs is based on the counting of relevant canonical dimensions. Dynamical 
models have two scales, the time scale T and the length scale L, consequently 
the physical dimension of any quantity F can be defined as [F] = L~ dp T~ dF , in 
which dp and d F are the momentum and frequency dimensions of F, In diffusion 
models, these dimensions are always related to each other since dt ~ d\ in the 
diffusion equation that allows us to introduce a combined canonical dimension, 
dp = d F + 2d F . One can check out that each term in (| 1 T[) is dimensionlcss 
if the following relations hold: dp = 0,d u = d, d v ,\ = —2 + 2-1 = 0, and 
d g = 2 — (d+ 1) . The field theory (fT7|) is logarithmic (the conventional dimension 
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Figure 4: The first diagrams for the inverse exact responde function of the magnetic 
field in MHD. The first diagram corresponds to the responde function of the linearized 
problem, Rij b , in l|28j) . The second diagram gives the second order (in the coupling 
constants) one-point correction due to fluctuations of the magnetic field. The third 
diagram stands for the first order two-points contribution. The forth diagram describes 
the the combined effect of the convection and the stretching. The fifth diagram gives 
the second order two-points correction due to the combined effect of hydrodynamic 
drag, convection, and stretching. 

of the coupling constant d 9i — 0) for the Burgers equation (d = 1), while in two 
dimensions d gi = — 1, and d 9i = — 2 for d = 3 (the NS equation). Thus, in the 
infrared region (small moments, large scales) the diagram series in <?i define just 
the corrections to the diffusion kernel as d > 2. However, in the case of Burgers 
equation, all diagrams look equally essential in large scales. 

The diagrams diverge in the ultraviolet region (large moments, small scales) 
if their canonical dimension 

d r = -d^N^ - dpNp > 

where and Np are the numbers of corresponding external legs in the graph 
r. For the Green function (|16[) . we have Np = and = 1. Therefore, there 
is no ultraviolet divergent graphs in the diagram series ([24]) . At first glance, it 
seems that any graph having no external <f>— legs = 0) and any number of 
auxiliary fields 0' [Np > 0) should diverge since dp = 0. However, such a graph 
is also convergent in small scales because of the derivatives in the vertex factors 
V(k) which are always taken outside the graph onto the external <f>'— legs that 
effectively reduces its canonical dimension to d' r = dr — Np < 0. Therefore, the 
field theory with the action functional (jTTJ) has no ultraviolet divergences and 
does not need a renormalization. 

Each Feynman diagram in corresponds to a certain magnetohydrody- 
namic process. For instance, the very first diagrams displayed in Figs. 13141 
present the solutions of diffusion equations (in d-dimensional space). They de- 
scribe the simple viscous dissipation of a hydrodynamic vortex with no bifurca- 
tions and the motion of magnetic field through the fluid, following a diffusion 
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law with the resistivity of the plasma A serving as a diffusion constant. 

The second graph (#2) in Fig. [3] corresponds to a bifurcation characterized 
by the twofold splitting of the moment, k = q + (k — q). Under the spatial 
Fourier transformation, it is equivalent to the following analytic expression: 

oc 

r x (k,t) = -gfdt'k- A(M-i') 

o (oU) 
/ (q, f ) v (q) A (k - q, t') v Q (k - q) , 

where A (k, t) is the spatial Fourier transform of the diffusion kernel and Vo (k) 
is the Fourier spectrum of initial condition. It is worth to mention that the 
diagram expansion for the Green function can be also discussed for the function 
defined on a finite domain supplied with periodic boundary conditions. In the 
latter case, it has a discrete set of harmonics, and the integral (j3"0)) turns into 
sums. 

The time integration in (|30[) can be performed easily, 



Tj — —g k • exp (- 



J, 



k><? 2 (27r) d 



vk 2 (t - t )) 

v (q) v (k - q) 2y(q \ k _^ ) 



(31) 



The singularities in T\ appear at q = and q = k, when the vortex does not 
bifurcate. The remaining momentum integral in (|31[) can be interpreted as an 
expectation value, 

oo 

Pi (k) = fdrf [v (q) v (k - q) exp {-2v (q ■ k - q 2 ) r)] , (32) 
J J (2ir) 
to 

over the Poisson process of vortex bifurcation at momentum k. Bifurcation 
of vortexes is the Poisson stochastic process developing with time [2]. Further 
diagrams shown on Figs. 13141 represent more complex magnetohydrodynamic 
processes. 



6 The multiplicative factors of Feynman diagrams 
in MHD 

One of fascinating features of the proposed stochastic approach to MHD is 
that all physically admissible behaviors of the MHD system are encoded by the 
certain integer solutions of the Eq. (f2"5]) . It is clear that not all such the solutions 
are equally contribute to the diagrammatic series . While drawing diagrams 
admissible with respect to the " grammar" prescribed by the standard Feynman 
rules discussed in the previous section, one can see that only some particular 
combinations of coupling constants can appear as the factors before diagrams 
for the Green functions. 
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In order to find these factors, we use the standard methods of graph theory. 
Indeed, each diagram constitutes a graph, in which every node representing one 
of four physical interactions is proportional to the relevant coupling constant. 
We start drawing a diagram from the point (x, f) by adding vertices and con- 
necting them in a way admissible with respect to the Feynman rules. At any 
step, the tails representing the fields v and b can be fused together at the ver- 
tices of the ultra- local interactions (the composite operators) v' (xo,io) and 
b' (x' Q ,t' ), m > 1. It is convenient to count the powers of coupling constants 
in the multiplicative factors before diagrams of perturbation theory with the 
use of the "grammar matrix" G given in the Appendix [B] 

The grammar matrix is the weighted connectivity matrix, which expresses 
the fact that the tails <f>, </>' belonging to the interaction vertices V^™, V^ a , 
and (being proportional to g\, 52, 53, and 54) can be connected to each 

other accordingly to the Feynman rules only by means of the propagators (|27p . 
For example, the first row of the grammar matrix G shows that the v'— tail 
belonging to the vertex V^ a v may be connected within a Feynman graph either 
to the v— tails in the similar interaction vertex V t ^" s v or to the v— tails in the 
interaction vertex V-^ . In both cases, the diagram amplitude acquires the fac- 
tor g%. It is then obvious that all possible diagram structures admissible by the 
Feynman rules can be reproduced by the powers of the grammar matrix, G ra , 
n > 1. Here n is the number of interaction vertices appearing in the diagram 
(excepting the ultra-local vertices v' (xo,t Q ) and b' (x' ,t' ) and the starting 
vertex related to the point (x, t)); it can be interpreted as the number of bifurca- 
tions of moments in the certain magnetohydrodynamic process. The entries of 
G n are the monomials, the products of powers of the coupling constants 51,2,3,4- 

The use of arguments given in Appendix[B]helps to verify that the multiplica- 
tive factors for Feynman diagrams of the Green functions G v and Gb contain 
powers of Gv = 9\929z9a, Gb = 925354 (see Tab. 1). Namely these combinations 
of coupling constants play the role of the expansion parameters in the diagram 
series (j2~4"|) for G v and Gb- 

Table 1: The multiplicative factors of Feynman diagrams in MHD 



# of bifurcations 


G v 


Gb 


n = 2k, k > 


9i-G k v , 
52 • Gv, 


53 -GZ, 
9 a -Gt, 


n = 2k+ 1, k> 


515254 • Gv, 
515253 • Gv, 
9 2 i92-G k v , 


g'iai-Gb, 

5I53 -Gt, 

52535 '4 • Gb, 



Consequently, if in the equation k\ + ki + k% + k^ — N we assume that 

1. either ki = 2fc+ 1, k% = &3 = £4 = k, or k\ = 2k, ki = k+ 1, k% = k^ = k, 
then the diagram series ([24]) reproduces the Feynman graphs for the Green 
function of fluid velocity G v containing an even number of bifurcation of 
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moments; 



2. either ki = 2k+l, fc 2 = &4 = k+1, k?, = k, or ki = 2fc+l, /c 2 = k% = k+1, 
ki = k, or eventually k\ = 2k + 2, fc 2 = k + 1, &3 = £4 = fc, then the 
diagram series (|24p reproduces the Feynman graphs for the Green function 
G v containing an odd number of bifurcation of moments. 

Diagrams for the Green function of magnetic field Gb do not contain the vertex 
responsible for the hydrodynamical interaction V^ s v , and therefore ki = 0. 
Then, if we suppose that 

1. either /c 2 — k, k$ = 2k + 1, fc 4 = 2k, or fc 2 = k, k% = 2k, k± = 2k + 1, then 
the diagram series (f24l) reproduces the Feynman graphs for Gb containing 
an even number of bifurcation of moments; 

2. either fc 2 = k, k 3 = 2k + 2, k t = 2k + l, or k 2 = k, k 3 = 2k + l, fc 4 = 2k + 2, 
or eventually k% = k + 1, k$ = 2k + 1, fc 4 = 2k + 1, then the diagram series 
(f2"4")l reproduces the Feynman graphs for the Green function Gb containing 
an odd number of bifurcation of moments. 

We have suggested everywhere that k E Z + . 

7 The large order asymptotic behavior for the 
MHD Green functions. Instanton approach 
and Borel summation 

It is important to note that we do not know apriori whether the coupling 
constants 31,2,3,4 are small or large. To get the information on the convergence 
of asymptotic series we should study the asymptotic behavior of large 

order coefficients Q^ 1 ^ 2 ^ 3 '^ ) j {1 _|_ fc 2 -f /j 3 + fc 4 = TV, in the diagram series as 
N — > 00, by the asymptotic calculation of the Cauchy integral, 

y~,(ki,k2,k 3 ,ki) __i 1 r tt-4 dg^ 



in which G^ (x, t; xo, to; x' , t' ) is the functional integral (|16p . The contour of in- 
tegration in the multiple integral ([33]) embraces the points {gi = 0}, i = 1, . . .4, 
in the complex plane. 

We estimate the functional integral ([33)) by the steepest descent method sup- 
posing that Yli=i ki = N is large (the instanton method). In so far, the instan- 
ton approach has been applied to various problems of stochastic dynamics, see 
[39]-[43]. In contrast to all previous studies, in the MHD system we do not have 
one coupling constant, but four. In the previous section, we have demonstrated 
that by applying the additional conditions for the integers ki to the equation 
Si=i ki = N in the series (|24[) . we can derive the diagram expansions for the 




x G^(x,t;x ,t ;x ,t' ) exp(-X\ =1 h log^ 



(33) 
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different Green functions being in the diverse statistical regimes (characterized 
by the even and odd number of momentum bifurcations respectively). 

Following the traditional instanton analysis, we perform the uniform rescal- 
ing of variables in the action functional (|17[) in order to extract their dependence 
upon N, 



9i -»• 9i/N, v -> v N, A -> A N, x -> x VN. ^ 

This keeps the action functional (JT7J) unchanged, thus each term acquires the 
multiplier N and then formally gets the same order, as the log gi in (|33p indepen- 
dently of the type of interaction i = 1, . . .4. The corresponding Jacobians from 
the numerator and the denominator of (|16p cancel. The saddle point equations 
are 

Vi + gi v (vd) vi — giv (bd^j b l ~v Av, = S (x — Xq) 8 (t - t ) , 
53 A (vd) hi — g t A (bd^j v~i — A Abi — S (^x — 5 (t — i, 
' + gi v (vd) v- - 33 A (b'dj bi + v Av! t = 



b 



K + 92 v (v'd) h + g 3 A (b'd) Vi + X Ab- = 0, 
vv t (vd) v- = g^ 1 , vv[ (bd^j h = g? 1 , 
Xbl(vd)bi=g^, Xb;(bd) 



(35) 



Vi = 9i l - 

The first two equations in ([35)) recover the original Cauchy problem for the 
MHD equations. The next two equations (#3 and #4) occur within the frame- 
work of the stochastic approach since we describe the microscopic dynamics 
in the MHD system as Brownian motion. The equations for the auxiliary 
fields are characterized by the negative viscosity and resistivity, and therefore 
(t > to) — b' (t > t'^j = 0. The last four equations in ([33)1 determine the 

saddle-point values of the coupling constnts that allows to exclude the interac- 
tion terms from the previous saddle-point equations and reduce the system (|35[) 

to 

v ' K v v = v ' 5 (x — xq) 5 (t — to) , 
h' K b b = b 7 6 ( x - x') 5(t-t' 



v 



vify = i, 

hK*b' = 1, 



(36) 



in which we have introduced the diffusion kernels, K v = —iuj + vp 2 and K b — 
—ioj + Xp 2 (in the (ui,p) Fourier space) and K*, K£ are their Hermit conjugated 
forms. 

Bifurcations of vortexes arisen due to the nonlinear interactions in the MHD 
system do not conclude into a critical regime, and therefore the time spectrum 
in the nonlinear model is the same as for the free diffusion equations, T cx 
L 2 , that is the reason for the branching processes are Poisson distributed with 
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the characteristic times 1/ {y k 2 ) and 1/ (Afc 2 ). One can see that the saddle- 
point configurations |v',b', v, b| which satisfy should be independent of 
the Poisson branching processes, and therefore, the solutions could exist before 
bifurcations start that is as t — > min (t , t^j . With the use of the power model 
for the Dirac delta function, 

S (x) = lim 



oir (x 2 + e 2 ) ' 



one can find that in the limit t — ► min^oiiojj Q36j) is satisfied by the fol- 
lowing radially symmetric solutions {4> = {v,b}, <ft = { v }) X = 

<l>(T,t) = y/(r-T ) 2 + (t-t ) 2 , <t>'(r,t) = H^<P(r,t), 

^r,*) = ^arctan(^), ^(r,t) = ^ ^^i^ . 

(37) 

In the first solution (|37|) , the auxiliary fields i/ and b' have poles as either t — > to 
or £ — * t i and then it follows from |35|) that 5* = 52 — an d .93 =54 =0, the 
point which definitely lays inside the integration contour in (|33[) . In contrast 
to it, in the last equation in (f37|) . v — > as t — > to an d & — > as < — > f 
respectively, and therefore 5* 2 3 4 00 that is definitely outside the integration 
contour. There is a subtle point in Eq. (|37p concerning the second solution 
since v(r,t = t ) = b(r,t = t' ) = H(0)/2, and the position of the saddle- 
point configuration charges 5*234 depends upon the conventional value for the 
Heaviside function of zero argument, H(0). While estimating the functional 
Jacobian, we had assumed following the standard convention [7] that H(0) = 0. 
Then, one can easily verify that in this case we also have 5* 2 3 4 — ► 00, so that 
being interested in the large order asymptotic behavior of the Green functions 
in MHD we do not need to take the second and the third solutions into account. 
Even if one takes H(0) ^ 0, and then < g\ 2 3 4 < 00, it is always possible 
to deform the integration contour in ([33)) in such a way to avoid g\ 2 3 4 to be 
encircled. Therefore, the first solution in (|37|) is the only one we need, and 
substituting it into (|35|) . we can obtain the microscopic power models for the 
coupling constants, 

t&^gSb 81 *-^' (38) 

Fields (f>, (j/, and the coupling constants 51,2,3,4 fluctuate around their saddle- 
point values (/>*, </>'*, and g* 2 3 4- By means of the standard shift of variables, 
<p = 0* +50, S(f>' = 4>'x+6(j)', 31,2,3,4 = <?* 2,3 4+^51,2,3,4; one makes them fluctuate 
around zero, so that S<fi(oo) = 0, Sip' (00) = 0. Moreover, if we assume that the 
MHD system is isotropic (i.e., there is no the global bias of the magnetic fields 
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and the conductive fluid is isotropic), then tr 6<f> = tr 5(f)' = 0, and therefore all 
fluctuations posses the central symmetry, 5(f) = 5(f>(r,t), 5(f)' — 8(j>'{r,t), the 
same as the saddle point configuration. 

The contours of integrations over the variables 531,2,3,4 now passes through 
the origin, and are directed there oppositely to the imaginary axis. The in- 
tegrals over 531,2,3,4 are conducted now on the rectilinear contours in com- 
plex planes (ioo,— ioo) (in accordance to the standard transformation of con- 
tours in the steepest descent method). At the turn of the integration contours 
5ffi,2,3,4 - * —i 531,2,3,4, where the multiplier (—i) appears so that the result 
q^mMM) j g a j wa y S rea j The contributions to the Cauchy integral (|33[) comes 
from the poles £31,2,3,4 = —g* 234 an d tends to zero as t — > min (to, t' ) . The 
values of the functional integrals on the saddle point configurations are pro- 
portional to the entire volume of the functional integration, they cancel in the 
numerator and denominator simultaneously. While calculating the fluctuation 
integral (|16p . we take into account that the first order contributions in 5(f> and 
5(f>' are absent because of the saddle-point condition. We also neglect the high- 
order interactions between fluctuations, 0(5</> 3 ), 0(5(p 4 ), etc. to arrive at the 
Gaussian functional integrals for the Green functions G v and Gb, 

G v '(r,5t) ~ S t->o [Gv) ax 

x JJ V 5(f) V 5(f)' exp -f tr [5(f)' + % 5(f) (±d r ) 5(f)] , 

G\ >(r,5t) -st^o [G b ) 2 27r 

x // V 5(f) V 5(f)' exp -f tr [5(f)' K ^5 (f) + ± 5<j> (± d r ) 6<p] , 

(39) 

in which Q v and Qb are the the expansion parameters in the diagram series l|24p 
for G v and Gb introduced in Sec. [6] in concern with Tab.l. Performing the usual 
rescaling of fluctuation fields 

5(f) -> 5(f>/ 'VN, 5(f>' -> 5(f)' /VN, 

we compute the Gaussian integral, with respect to 5(f) first, and then the resulting 
Gaussian integral over the fluctuation of the auxiliary fields 5(f)' , 

G v N \r,5t) N N -V* exp(-f \og{G*)) det^ 1 (l + 0(£)) , 

G ( b N \r,6t) -^o^" 1 / 2 exp(-f log (^)) det^ 1 (l + 0(±)). 

(40) 

The kernels of the operators 1 are the Green function of the linear diffusion 

equations. Using the Stirling's formula, one can check that the coefficients G^ 
of the asymptotic series (|24|) demonstrate the factorial growth (like in the most 
of quantum field theory models): 

G v N \r,5t) ^expiV(l-i log (^)) det if^ 1 (l + O (^)) , 

G ( b N) (r,5t) ^expJV(l-|log(S 6 *))detJf^(l + 0(^j). 

(41) 
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Therefore, the asymptotic series can be summed by means of Borel's pro- 
cedure. It consists of the following transformation of series ([24]) 



EivGf^ = E,r(« + i)G^g; 

= EiV Jo dT G ( T ^) e > 

where G,f ] = G$° /T(N + 1) are the new expansion coefficients which do not 
exhibit the factorial growth. 

It is traditional, while performing the Borel summation, to change the orders 
of summation and integration in (|42p 

roo 

" N (43) 

N " u " ; A 

The we sum over N in the r.h.s. of (|4"3"1) . 

1 — 1 /.oo oo \iV 



e / ^ ^ (^) N e- T = dr e ~ r e «w 

« «/o Jo at 



17 " jV=l 



and obtain 



det iTT 1 



, dr e~ T log f 1 - t§| | . (44) 

2^ Jo \ s;J K ' 

The integration of (|4"4"1) over r gives us 

G< ^ ode tV(l + ^i(|)exp(-|)) (45) 

where Ei(ir) is the exponential integral defined as Ei(x) = — y^ 1 e~ v dy 
understood in terms of the Cauchy principal value at y = 0. Then we can 
use the microscopic models (|3"51) and recall that by definition 31,2 = 1/f and 
§3,4 = 1/A in order to estimate the both ratios Gt/Gv and Gl/Gb in P5| as 

G*a St 5 / 2 

t =*" V- < 46) 

In such a simplified model, the non-Maxwellian corrections to the distribution 
functions arisen due to the kinetic effects are accounted by the new distribution 
function. 



G^(6t,6r) ^st-*o 



ex p(-i|fe)) 

(47T X (^)) d/2 



i + ±.m exp -M^ 



(47) 
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8 Discussion and Conclusion 



Speaking rigorously, the classical MHD equations cannot be applied if the 
time scale of collisions is comparable or longer than the other characteristic 
times in plasmas. The particle distributions are far from being of the Maxwellian 
shape and the certain kinetic models giving an insight into the collision statistics 
have to be taken into account. The simplest such model which allows accounting 
the kinetic effects being nevertheless completely in the framework of the classical 
approach based on the MHD equations is suggested in the present paper. 

We investigate a simple generalization of the MHD model ([T|) modelling 
fluctuations of the configurations {v, b} considered as the stochastic fields (or 
the trial trajectories of the MHD system) for which the classical MHD (TT|) plays 
the role of the mean field equations. 

The essential point of our approach is that we have used the field theory 
formulation of the dynamics which allowed us the implementation of various 
powerful technics borrowed from the quantum field theory. In particular, with 
the use of the instanton technique and Borel's summation, we have computed the 
asymptotic series for the Green functions accounting for the kinetic effects as the 
corrections to the unperturbed diffusion kernel describing the pure relaxation 
dynamics. Similarly to the most of quantum field theory models, the high order 
contributions into the Green functions exhibit a factorial growth. 

It is interesting to compare the perturbed diffusion kernel (|47l) in the 3D 
space, d = 3, with the standard diffusion kernel depicted by a Gaussian curve. 
In Figs. 18181 we have sketched the profiles of standard diffusion kernel (the solid 
lines) calculated at at several consequent time steps for v = 0.2 together with 
the perturbed kernel profiles given by (j4"T|) . 

It is clearly seen that the essential corrections to the standard diffusion kernel 
are arisen in short times and small scales, while they are negligible in long times 
large scales. 
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Figure 5: The profiles of standard diffusion kernel G(8r) (the solid line) calcu- 
lated at at several consequent time steps t > for v = 0.2. The dash-dot lines 
present the asymptotic kernel (|4T[) (as t — » to) accounting for the perturbation 
due to the kinetic effects modelled by Brownian motion. 
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Figure 6: The profiles of standard diffusion kernel G(5t) (the solid line) cal- 
culated at several distant points for v = 0.2. The dash-dot lines present the 
asymptotic kernel (I47[) (as t — > to) accounting for the perturbation due to the 
kinetic effects modelled by Brownian motion. 
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A Analytical expressions for Feynman diagrams 

Below, we write down the analytical expressions (in ti-dimensional Fourier space 
(k, t)) correspondent to diagrams shown in Figs. (j3!4|) . 



G„(k, t - t ; k', t-t' ) = v (k)R v ' v (k, t - t ) 

+ S™dt'V v ' m {k)Ry'*{k,t-t') J -^Rv'v^t' ~ t )v (q)R v ' v (k- q,t' - t Q ) 
w (k- q) 

+ f °° dt'V*'™(k)R«'*(k, t-t')f j^R v ' v (k - q,f - t )v (k - q) 

/ dt"V v '™(q)R"'*(q,t' -t") J J 0^R^(p, t '-t o )v o (p)R v ' v (q-p,t'-to) 

uo(q-p) 

+ f Q °° dt>V*' bb (k)R^(k, t-t')J j^R b ' b (q, t> - t' )b (q) 
R b ' b (k~q,t'-t' )b (k-q) 

+ / °° dt"V*' m (k)fl»'*(k, t" - t) / ^?«o(k - p)iT'"(k - p, t" - i ) 
jj~ dt / 7 »'66( p)ir ,'t, (pj f / _ n J ^R^q, t> - t' )b (p)b (p - q) 

R b '\p-q,t' -f Q ) 

+ J Q °° dt"V v ' vv (k)R v ' v (k, t - t") J ^Rv'^p, t' - t")R v ' v (k - p, t'" - t") 
V*' bb (p)V v '™(k - p) / j^R b ' b (q, ? - «&)6o(p)6o(p - q) / ^^(s, f" - i )) 
i?"'"(k p s, f' - to)«o(s)«d(k - p - s) 

+ / °° dt"V v ' m (k)R v ' v (k, t - t") f °° dt'" f j0^R b ' b {p, t'" - t")V b ' bv {p) 
V b ' b -(k - p)R b '»(k - p, t" - t>) J T ^v ( q )vo^q)bo(q)bo(-q)R b ' b (q, t'" - t' ) 
R?' v {q,t' - t Q )R v ' v (-q,t"i - t a )R b ' b (-q,t' - t' ) 

(48) 

G b (k, t-t ;k',t- t' Q ) = b (k')R b ' b (k', t - t' ) 

+ f Q °° ^'^'^(kO^'^k', t-t')J j^R v ' v (k' -q,t'- f )6o(k - q) 

I dt"V-' bb (q)R b ' b (q, if - t") J j^R^ip, t> - t' )bo(p)R b ' b (q - p,f - t' ) 

feo(q-p) 

+ / °° dt'V b ' vb (k')R b ' b (k', t-t')J T ^ J v (s)R"'%k' - s,f - t ) 
J^bo(q)R b ' b (k'-q,t'-t' ) 

+ J Q °° dt'V b ' vb (k')R b ' b (k\ t - t') f Q °° dt" J J 0^ bo (s)R v '-(k' - s, t" - t' ) 
J j^v (k' - s - q)v (k' + s)R v ' v (k' -s-q,t'- t Q )R v ' v (k' + s,t" - to) 
V b ' vb (k' q s)R b ' b (k' - s - q,t' - t") 
+ Jo°° dt " Io° dt'V b ' vb (k')R b ' b {k', t - f ) 

J j^V b '^(k' - q - s)R b ' b (k' - 8 - q, t< - t") J^J jgpR^iq, t> - t )v (q') 
v (k' - q')R v ' v (k' - q',t" - t )R v ' v (q', t" - t )V v ' vv {k! - p) 

(49) 
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B The "grammar" matrix for the Feynman di- 
agram technique 

The "grammar" matrix G is the weighted connectivity matrix, which expresses 
the fact that the tails 0, </>' belonging to the interaction vertices V^ a v , , 
and (being proportional to gx, g%, 53, and g±) can be connected to each 

other accordingly to the Feynmann rules only by means of the propagators (f2"T|) . 
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The starting vertex in a diagram contributing into the Green function G v 
may be either V^ m v s v or V^f. In the first case, the multiplicative factor acquires 
the additional multiplier gi, and it is gi if the second vertex is used as the 
starting one. In order to include the starting nodes into account, we introduce 
the vector T = [gi, 0, 0, 172: 0, 0, 0, 0, 0, 0, 0, 0]. Then, the multiplicative factors 
of the diagrams for the Green function G v are written as following, 

[Gy factor - {GM„>o- (50) 

Similarly, we can find the multiplicative factors arising in the diagrams for the 
Green function Gb by 

[^factor = {<&"*>}„>„, (51) 

in which the vector b T = [0, 0, 0, 0, 0, 0, #3, 0, 0, #4, 0, 0] expresses the fact that 
V-^s is the starting vertex for G&. 
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